Isoscaling of the Fission Fragments with Langevin Equation^ 
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Langevin equation is used to simulate the fission process of 112 Sn + 112 Sn and 116 Sn + 116 Sn. The 
mass distribution of the fission fragments are given by assuming the process of symmetric fission or 
asymmetric fission with the Gaussian probability sampling. Isoscaling behavior has been observed 
from the analysis of fission fragments of both reactions and the isoscaling parameter a seems to be 
sensitive to the width of fission probability and the beam energy. 

PACS numbers: 24.75.+i, 25.85.Ge, 21.10.Tg 



Since the isoscaling law has been observed experimen- 
tally 0, 0, Q , many statistical models have successfully 
explained the isoscaling behavior. Isoscaling means that 
the ratio of isotope yields from two different reactions, 
1 and 2, R 21 {N,Z) = Y 2 {N,Z)/Y X (N,Z), is found to 
exhibit an exponential relationship as a function of the 
neutron number N and proton number Z Q 



R 21 (N, Z) 



Y 2 (N, Z) 
Yi{N,Z) 



Cexp(aN + (3Z). (1) 



where C, a and (3 are three parameters. In grand- 
canonical limit, a — A[i n /T and f3 = A/x z /T where 
A[x n and A[i z are the differences between the neutron 
and proton chemical potentials for the two reactions, re- 
spectively. This behavior is attributed to the difference 
of two reaction systems with different isospin asymme- 
try. It is potential to probe the isospin dependent nu- 
clear equation of state by the studies of isoscaling Q. 
So far, the isoscaling behavior has been experimentally 
explored by various reaction mechanisms, ranging from 
the evaporation 0] , fission |f| and deep inelastic re- 
action at low energies to the projectile fragmentation 
and multi-fragmentation at intermediate energy 
Hcj . While, the isoscaling phenomenon has been 
extensively examined in different theoretical frameworks, 
ranging from dynamical model, such as BUU model 
and anti-symmetrical molecular dynamics model [ill ], to 
statistical models, such as the expansion emission source 
model, statistical multi- frag mentation model and the lat- 
tice gas model H El HI El Hi- 

In this work, we present an analyse for the fragments 
from the fission which was simulated by Langevin equa- 
tion. The isotopic or isotonic ratios of the different frag- 
ment yields from 116 Sn + 116 Sn and 112 Sn + 112 Sn system 
are presented and the features of the isoscaling behavior 
in fission dynamics are investigated. 

The process of fission can be described in terms of col- 
lective motion using the transport theory [Isl Hfil ItR Hij l . 
The dynamics of the collective degrees of freedom is 
typically described using the Langevin or Fokker-Planck 



equation. In this Letter, we deal with a Combine Dy- 
namical and Statistical Model (CDSM) which is a com- 
bination of a dynamical Langevin equation and a statis- 
tical model to describe the fission process of heavy ion 
reaction. This model is an overdamped Langevin equa- 
tion coupled with a Monte Carlo procedure allowing for 
the discrete emission of light particles. It switches over to 
statistical model when the dynamical description reaches 
a quasi-stationary regime. We first specify the entrance 
channel through which a compound nucleus is formed, ie. 
the target and projectile is complete fusion. 

In this work the total initial excitation energy E% ot is 
given by E* tot = Ei ab A T /(A T +A P )+Q where A T and A P 
represents the mass of target and projectile, respectively, 
and Q is the fusion Q-value calculated by Q = Mt + 
Alp — Mq^. Mt and Mp is the mass of projectile and 
target come from experimental data, respectively. If it is 
unavailable, it is calculated by macroscopic-microscopic 
model . Mq® is the mass of the compound nucleus 
which is calculated from the liquid-drop model. 

The dynamical part of CDSM model is described by 
Langevin equation which is driven by the free e nerg y F. 
F is related to the level density parameter a(q) 



F(q,T)=V(q)-a(q)T 2 



(2) 



in the Fermi gas model, where V(q) is fission potential. 
The overdamped Langevin equation reads 



dq 
~dt 



1 



dF{q,T) T 



Mf3(q) y Oq 



) + y/D{q)T(t), (3) 



where q is the dimensionless fission coordinate defined 
as half of the distance between the centers of masses of 
the future fission fragments. T(t) is a time-dependent 
stochastic variable with Gaussian distribution. Its aver- 
age and correlation function is written as 



< r(t) >=0, 
< T(t)T(t') >= 2S £ (t - t'). 



(4) 



The fluctuation strength coefficient D(q) can be ex- 
pressed according to the fluctuation-dissipation theorem: 



2 



D(q) = 



T 



Mf3 (q)> 



(5) 



where M is the total mass and 0o(q) is the reduced fric- 
tion parameter which is the only parameter of this model. 

The potential energy V(A, Z,L,q) is obtained from the 
finite-range liquid drop model |21j 



V(A, Z, L, q) = a 2 [l - k( 



N 



A 



2 ]A m [B s {q)-l] 



+c 3jm [B c ( q ) - 1] + c r L 2 A- 5 / 3 B r (q), (6) 

where B s (q), B c (q) and B r (q) means surface, Coulomb 
and rotational energy terms, respectively, which depends 
on the deformation coordinate q. a 2 , C3, k and c r are 
parameters not related to q. In our calculation we take 
them according to Ref. flif . 

We use c and h H3| to describe the shape of nucleus, 



1 



b 



p\z) = (1 - -.)((- - -Zy + B sh (c,h)z% (7) 



where 



Co 



cR, i? = I.I6A 1 / 3 . 



(8) 



The nuclear shape function B s h (c, h) and the collective 
fission coordinate q(c, h) of mass number A is expressed 
as 



B sh (c,h) = 2h- 



c- 1 



q(c,h) = |c(l + ^rS s h(c, /i)c 3 ) 



(9) 



The fission process of Langevin equation is propagated 
using an interpretation of Smoluchowski |23l| . In our cal- 
culation we adopt one-body dissipation (OBD) friction 
form factor /?o((?) [3 which is calculated with one-body 
dissipation with a reduction of wall term. Here we use an 
analytical fit formula which was developed in Ref. [25| 



f3oBD(q) = { 



15/q' 



0.43 



1 - 10.5<7 u ' 9 + <f if q > 0.38 



32 - 32.21<7 



if q < 0.38 



In the dynamical part of the model the emission of light 
particles (n,a,p,d) and giant dipole 7 are calculated at 
each Langevin time step r, the widths for particle and 
giant dipole 7 decay are given by the parametrization of 
Blann [2(| and Lynn j2?|], respectively. 

Within the framework of Langevin simulation we chose 
200,000 fission events which happen on dynamic chan- 
nel (we give up the events which happen in statistic 
part of CDSM model) and chose a Gaussian distribu- 
tion random number as the mass asymmetry parameter 
«o = (Ai — A2)/(Ai + A2), which is defined as the ra- 
tio of the volumes of two parts of the nucleus obtained 



when it reaches the scission point. When ag = means 
symmetrical fission. It is taken from a Gaussian distri- 
bution random number from -1 to 1 with the mean is 0. 
Ai and A2 is the mass of the two fission fragments, re- 
spectively. In this work we assume the fission fragments 
have the same N/Z as the initial system and then Z\ or 
Z2 of fission fragments can be deduced from A\ or A2. 
This assumption is similar to the case of of deep inelastic 
heavy ion collisions at low energies, where the isospin de- 
gre e of freedom has been found to first reach equilibrium 

From a practical point of view, the isoscaling occurs 
when two mass distributions for a given Z from two pro- 
cesses with different isospin are Gaussian distributions 
with the same width but different mean mass. Basically, 
the isotopic distribution can be described by 



Y(N,Z)=f(Z)exp[ 



(N - N z f 
24 



(10) 



where Nz is the centroid of the distribution and <t 2 z de- 
scribes the variance of the distribution for each element 
of charge Z. This leads to an exponential behavior of the 
ratio R21 if the quadratic term in Nz is neglected, 



ln{R 21 ) 



[{Nzh ~ (N z )i]N 



(11) 



Note that Eq. l|ll|) requires the values for the cr 2 z to be 
approximately the same for both reactions, which is a 
necessary condition for isoscaling. Indeed, we observed 
this case in our simulations for Sn + Sn collisions. In the 
Langevin equation, <j 2 z mainly depends on the physical 
conditions reached, such as the temperature, the density 
and the friction parameter etc. Considering that R' 21 — 
exp(aN) for a given Z, a ~ [(jVz)2 ~ (jVz)l1 . Usually oi 
can be considered to be proportional to temperature T 
of the fragments, in this way 



[{Nzh ~ ( N z)i 
T 



(12) 



where [{Nz)2 — (^z)i] can be understood as a term of 
the chemical potential difference between two reactions. 

Eq.(l) can be written as Zni? 2 i — Cz + o>N, where 
Cz = InC + fiZ, if we plot R21 as function of N, on a 
natural logarithmic plot, the ratio follows along a straight 
line. In Fig. ^ this behavior is observed in Langevin sim- 
ulation. From the figure, the relationship between a and 
the charge number Z of the fission fragments can be de- 
duced. In order to investigate the effect of the width 
of Gaussian probability distribution on the isoscaling be- 
havior, we change the widths of the Gaussian distribution 
of mass asymmetry parameter a>o, namely a ao = 0.04, 
0.06, 0.08 and 0.20, with the random number from -1 to 
1 and the mean of 0. Fig[21 shows the isoscaling param- 
eter a as a function of Z in the conditions of different 



3 



Cq - From this figure, we know in the low c Qo , i.e., the 
fission fragments are overwhelmingly dominated by the 
symmetric fission, a increases with Z. This means that 
the isospin effect becomes stronger with the increasing 
of Z. In a recent analyse of Friedman [f| with a simple 
liquid-drop model, a systematic increase of the isoscaling 
parameter a with the proton number of the fragment ele- 
ment is predicted. In our simulation, this kind of increase 
of a with Z stems from the dominated symmetric fission 
mechanism. While, in the another extreme case from the 
Fig|21 fe. with the larger a ao , a shows a contrary trend 
with Z, i.e., it drops with Z. In this case, it seems that 
there exists stronger isospin effect for the fragments with 
lower Z. In a middle case, the rising branch and falling 
branch competes with each other, the mediate isoscaling 
behavior appears and a minimum of a parameter occurs 
around the symmetric fission point. We note that the 
fission data of 238 > 233 u targets induced by 14 MeV neu- 
trons reveals the backbending behavior of the isoscal- 
ing parameter a around the symmetric fission point Q 
as stated above. In their study, they interpreted this 
originates from the temperature difference of fission frag- 
ments since the isoscaling parameter is typically, within 
the grand-canonical approximation, considered inversely 
proportional to the temperature ( a = A[i n /T ) as stated 
above. It is, however, not a priori obvious why such a 
grand-canonical formula can be applied to fission. In our 
case, this kind of backbending of isoscaling parameter 
apparently stems from the moderate width of the fission 
probability of the fissioning nucleus as shown in Fig|2 
Of course, we did not exclude the change of temperature 
of the fission fragments due to the change of the width of 
the fission probability. Overall speaking, we find that the 
isoscaling behavior is sensitive to the width of the fission 
probability distribution. 

In addition, the simulations are systematically done in 
different beam energies. The values of a are extracted as 
a function of beam energy for the fragments Z = 44-52 as 
shown in FigEl It shows that a decreases as the beam en- 
ergy increasing which means that the isospin effect fades 
away with the increasing of Ei a b- This behavior is similar 
to that a drops with the temperature in the statistical 
models as well as experiments 0, 0, HjJ • 

In summary, we applied Langevin model to investigate 
the isoscaling behavior in the dynamical process of com- 
pound nuclear fission. In order to treat the fission frag- 
ments, we assume that the mass asymmetry parameter 
of the fissioning nucleus is taken from a random number 
of Gaussian distribution whose width is o~ ao . The simula- 
tion illustrates that the yield ratios of fission fragments in 
the dynamical fission of 116 Sn+ 116 Sn and 112 Sn+ 112 Sn 
reaction system shows the isoscaling behavior. More in- 
terestingly, we found that the isoscaling parameter a is 
sensitive strongly to the Gaussian width tr Qo of the mass 
asymmetry parameter. When o~ ao is small, i.e. the fis- 
sion is almost symmetrical, a increases with the atomic 




FIG. 1: The yield ratio of the fission fragments between 116 Sn+ 116 
Sn and 112 SVi+ 2 Sn in the Langevin model with a aQ = 0.06 and 
E/A = 8.4 MeV. Different symbols from left to right represent the 
calculated results for the isotopes from Z = 37 to 59. The lines 
represent exponential fits to guide the eye. 




FIG. 2: The scaling parameter a as a function of Z in the different 
width (<T ao ) of the mass asymmetry parameter cto with Gaussian 
distribution for fission. 



number of fission fragments, which is similar to the the- 
oretical prediction of a simple liquid-drop model |5J. In 
contrary, when a aa is large, for instance, a ao — 0.20, a 
drops with Z of fission fragments. However, in the inter- 
mediate values of a ao , a shows a backbending with Z of 
fission fragments, which is similar to the observation of 
the 238 < 233 TJ fission data induced by 14 MeV neutrons [||. 
In addition, it is found that a drops with the beam energy 
of the projectile, reflecting the temperature dependence 
of isoscaling parameter. In general, the isoscaling anal- 
ysis of the fission data appears to be a sensitive tool to 
investigate the fission dynamics. 
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FIG. 3: a as a function of beam energy for the fragments Z 
44-52. The width of the Gaussian probability a ao is 0.06. 
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